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Abstract 

The static force Fsir) and the strong couphng apir), which defines the gluon- 
exchange part of Fsir), are studied in QCD background perturbation theory 
(BPT). In the region r < 0.6 fm ai?(r) turns out to be essentiahy smaller than 
the coupling aB(r) in the static potential. For the dimensionless function $_B(r) = 
r^Fe(r) the characteristic values ^b{i"i) = 1-0 and ^b{^o) = 1-65 are shown to be 
reached at the following QQ separations: rl^/a = 0.77, ro^/o^ = 1.09 in quenched 
approximation and r\yfa = 0.72, tq^/g = 1.04 for nj = 3. The numbers obtained 
appear to be by only 8% smaller than those calculated in lattice QCD while the 
values of the couplings apiri) and apiro) in BPT are by ~ 30%(nj = 3) and 
50% (nj = 0) larger than corresponding lattice couplings. With the use of the 
BPT potential good description of the bottomonium spectrum is obtained. 

1 Introduction 

The static QQ interaction plays a special role in hadron physics and 
the modern understanding of the spin-independent part of the static 
potential V{r) comes from different approaches: QCD phenomenology 
|T] - |3] , perturbative QCD IH [S] , the analytical perturbation theory [Bj , 
the background perturbation theory (BPT) [7|,|B|, and lattice QCD P- 
jTTj. It has been established that in the region r > 0.2 fm this potential 
can be taken as a sum of nonperturbative (NP) confining potential plus 
a gluon-exchange term: 

V{r) = VNp{r) + VGE{r). (1) 

Such an additive form was introduced by the Cornell group already 30 
years ago, just after discovery of charmonium [Ij, and later it has been 
supported by numerous calculations in QCD phenomenology and also on 



the lattice, where for the lattice potential the parametrization like the 
Cornell potential is used [9,10], 

VUr) = ar-^-^, (r > 0.2 fm) (2) 
r 

with eiat = const, aiat = \Giat- (For the lattice static potential we use 
here the notation aiat{e-iat) to distinguish it from the phenomenological 
case). 

In BPT the additive form ([TI) can automatically be obtained in the 
lowest approximation when both NP and Coulomb terms satisfy the 
Coulomb scaling law with an accuracy a few percents, in good agree- 
ment with the lattice data 

The explicit form of NP term is well known now: VNp{r) = ar at 
the QQ separations r > Tg, where Tg is the gluonic correlation length, 
Tg = 0.2 fm, measured on the lattice [13J. Such linear behavior is valid up 
to the distances r ~ 1.0 fm, while at larger r,due to the qq-paiv creation, 
the flattening of confining potential, is expected to take place [14]. Here 
in our paper we restrict our consideration to the region r ^ 1.0 fm. 

Numerous calculations of the meson spectra have shown that linear 
potential defines gross features of the spectra, in particular, the slope and 
the intercept of Regge - trajectories for light mesons [15]. At the same 
time the splittings between low-lying levels, the fine-structure splittings, 
and the wave function at the origin in heavy quarkonia are shown to 
be very sensitive to the gluon-exchange part of the static potential (jT]) 

V{r)=ar + VGE{r), (3) 

VcEir) = -3^— • 

Unfortunately, our knowledge of this term (or the vector coupling ast{r)) 
in IR region remains insufficient and the choice of ast{r) essentially dif- 
fers in different theoretical approaches. The only common feature of 
ast(r), used in QCD phenomenology, BPT, and supported by the lattice 
calculations, is that this coupling freezes at large distances, 

ast{ large r) = const = acnt- (4) 

but the true value of acrit is not fixed up to now. For example, in the 
Cornell potential, used in the phenomenology, the typical values of the 
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Coulomb constant in (I2j) lie in the range 0.39-^0.45 [1-3,16] while on the 
lattice in the same parametrization (Ej) aiat has the value by 50-^30% 
smaller: in quenched approximation (n/ = 0) aiat = 0.20 -!- 0.23 [^-[H] 
and in (2 + 1) QCD = 0.30 [H]. 

Even larger freezing value is obtained if the asymptotic freedom (AF) 
behavior of ast{r) (or ast{q) in momentum space) is taken into account. 
In BPT large acrit = 0.58 ± 0.02 follows from the analysis of the split- 
tings between low-lying levels in bottomonium [16,17] and this number 
is in striking agreement with acru, introduced by Godfrey, Isgur in [3] 
in phenomenological way. From the analysis of the hadronic decays of 
the r-lepton the number as{lGeV) = 0.9 ± 0.1 was determined in Ref. 
[18], while essentially smaller value, as{lGeV) = 0.45, was obtained 
in perturbative QCD with higher order corrections [19]. The Richard- 
son potential as well as "the analytical perturbation theory" give even 
larger acrit{nf = 3) = 1.4 [6j. Thus at present the true freezing value 
of the vector coupling, as well as ^^(l GeV), remains unknown. 

Besides the vector coupling ast{r) (in the static potential) the coupling 
apir), associated with the static force, can be introduced. 



Then an important information about the derivative ^^^(r) can be ob- 
tained from the study of the force. Note that the matrix elements over 
the static force define the squared wave function at the origin for the nS 
states in heavy quarkonia. In the lattice analysis the vector coupling aiat 
in Viat{r) (ED is supposed to be independent of the QQ separation r over 
the whole region 0.2 fm -^1.0 fm [10], i.e. a'l^^ = 0, and therefore in this 
region aiat = ap = const. This statement is in accord with recent calcu- 
lations of the lattice force in [11] where equal values a^*(ri) = a^*(ro) 
at the points ri ~ 0.34 fm, ro = 0.5 fm have been obtained. 

Different picture takes place in BPT where the vector coupling as(r) 
at the points ri = 0.34 fm and tq = 0.5 fm changes by 20% (see Section 
4) being close to the freezing value only at large distances, r > 0.6 fm. 




(5) 




(6) 



where by definition apir) is expressed through as<(r). 




(7) 
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Also due to the r-dependence an essential difference between ^^(r) and 
asir) occurs just in the range 0.2fm < r < 0.6 fm. 

Therefore several features of the vector coupling still need to be clar- 
ified. 

First, why in lattice QCD and in QCD phenomenology the Coulomb 
constant, used in the same Cornell potential, differs by 30% {uf = 3)? 

Second, whether the r-dependence of the vector coupling ast{r) (also 
of apir)) really exists and why it is not observed on the lattice? 

Third, what is the true freezing value? 

To answer some of these questions we shall use in our analysis here 
BPT which gives a consistent analytical description of the vector coupling 
both in momentum and in coordinate spaces. To test the BPT conception 
about the vector-coupling behavior in the IR region, recently the heavy- 
quarkonia spectra have been successfully described in this approach with 
the use of only fundamental quantities: the current (pole) quark mass, 
Kj^ijif)^ and the string tension [16], [17]. It is important that in BPT 
the vector coupling asiq) has correct perturbative limit at large g^, and 
therefore it is fully defined by the QCD constant A;^^ and also by so- 
called background mass Mb which is proportional y/a. 

Additional information about QQ static interaction can be extracted 
from the study of the static force Fb{t) in BPT, with further comparison 
to recent lattice results from [11]. To this end it is convenient to calculate 
the dimentionless function $(r) = r^F(r) at two characteristic points - 
r\\fo^ ri^^^ where 

$(riV^) = 1.0 and ^{r^yf^) = 1.65. (8) 

On the lattice the values ^/o'r^i\ y/arQ^ are shown to be different in 
quenched approximation and in (2 + 1) QCD. The same function ^^(r), 
calculated here in BPT, aquires the values (8) at the points ^/ari, ^/o^rQ 
which appear to be only by 8% smaller than those on the lattice. How- 
ever, the freezing value acrit in BPT is shown to be essentially larger 
than aiat- To compare the BPT and the lattice potentials we calculate 
here the splittings between low-lying states in bottomonium and demon- 
strate that the lattice static potential cannot provide good agreement 
with experiment, while in BPT such an agreement takes place for the 
conventional values of the pole mass of b quark, the string tension, and 
Aj^{nf = 5). We give also explanation why between the freezing values 
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in BPT and in the Cornell potential, used in phenomenology, there is the 
difference about 25%. 

The paper is organized as follows. In Section 2 we shortly present 
the vector coupling in BPT and discuss the correct choice of the QCD 
constant Ay in the Vector- scheme. In Section 3 the characteristics of the 
lattice force are presented while in Section 4 the values r\^J~o ^r^^ are 
calculated in BPT. The difference between the vector couplings in both 
approaches is also discussed. In Section 5 the splittings between low-lying 
states in bottomonium are used as a test to compare the static potentials 
from the lattice data, in phenomenology, and in BPT. We show that the 
lattice static potential cannot provide good agreement with experiment. 
In Section 6 we present our Conclusion. 



2 The strong coupling in BPT 

In BPT the gluon-exchange term V^(r) defines the vector (background) 
coupling aBif) in the same way, as "the exact coupling" OLstif) is defined 
inEq. (3): 

y^(r)=ar + yJi;(r); yj^^(r) = (9) 

With the use of the Fourier transform of the potential VqeW) back- 
ground coupling in coordinate space can be expressed through the cou- 
pling in momentum space [20]: 

2 /""^ dq 

asir) = - —sin{qr)aBiq), (10) 
7^ Jo Q 

where the vector coupling in momentum space is defined at all momenta 
in Euclidean space and has no singularity for q^ > [7]. In two-loop 
approximation the coupling 



g2 + 



t5 = ln^^-^, (12) 

contains the background mass Mb which enters under logarithm as a 
moderator of the IR behavior of the perturbative coupling. In asiq) 
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the Landau ghost pole disappears while the short distance perturbative 
behavior, as well as the Casimir scaling property of the static potential, 
stays intact [12]. 

The value of Mb is determined by the lowest excitation of a hybrid 
[21], however this mass cannot be considered as an additional (fitting) 
parameter in the theory, since in QCD string theory it can be calcu- 
lated on the same grounds as mesons, being expressed through the only 
(besides Ay) dimensional parameter ^/a•. 



We suppose here that in the static limit the coefficient ^ does not 
depend (or weakly depends) on the number of flavors nj. Direct cal- 
culation of Mb has not yet done and therefore the number ^ has been 
extracted from two fits: from the comparison of lattice static potential 
to that in BPT [20] and from the analysis of the spectra in charmonium 
and bottomonium [16,17], with the following result: 



In particular, the values Mb = 1.0 GeV and 0.95 GeV correspond to 
a = 0.20 GeV2 and 0.18 GeV^. 

Note that the logarithm ( IT^ in asiq) formally coincides with that 
suggested in Refs. [22] many years ago in the picture where the gluon 
acquires an effective mass nig and, as a result, in ( ITTl) instead of the 
value {"InigY enters. However, the physical gluon has no mass and in 
BPT the parameter Mb has been interpreted in correct way as a hybrid 
excitation of the QQ string, which is proportional y/a [21]. 

The background coupling aB{q) has correct PQCD hmit at ^ 
M^l and therefore the constant Ay(n/) (in Vector-scheme) under the 




logarithm ( |T^ can be expressed through the conventional QCD constant 
Aj^iuf) as in PQCD [23]: 



Mb = ^V^ 



(13) 



Mb = 2.236(ll)v^. 



(14) 






with ai = y — yn/, /5o = 11 — |n/. At present the values of are 
well established in two cases - from high-energy processes for n/ = 5 [24] 
and in quenched approximation from lattice calculations [25] : 




A^(2 - loop) = (216 ± 25) MeV , (16) 
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which corresponds to the "world average" as{Mz) = 0.117 ± 0.002, and 
the value 

A^(2 - loop) = ^-^^ = (237 ± 19) MeV (17) 

(with ro = 0.5 fm = 2.538 GeV"^) was calculated on the lattice in [25]. 

Then from ( ITBI) the corresponding values of Av{nf) in the Vector- 
scheme are following, 

= (295 ±35) MeV, (18) 

A^y^ = (379 ±30) MeV. (19) 

It is worthwhile to notice that the background coupling ^^(r) (12) is a 
universal function of the ratio 

AK) . M^. (20) 

and actually depends on the dimensionless variable x = rMg and A: 

2 f dx 

asir, Ay, Mb) ^ aBirMB', X) = - / — sm(x, rMB)aB(x, X). (21) 

7T J X 

The same ratio A also defines the freezing (critical) value of the coupHng 
(which coincides in momentum and coordinate spaces), and it is given 
by the analytical expression. 



47r ( Pilnt 



acritiq ^ 0) = acrit{r oo) = <! 1 - ^^^p^ \ , (22) 

terzi = ln-f = lnA2(n/). (23) 

Taking Ay from ( ITHl) . ( IT!1 and Mb (14) one obtains the following numbers 
for acrit- 

in quenched approximation 

aUnf = 0) = 0.419 (24) 

while for n/ = 5 the freezing value is by ~ 30% larger, 

+0 05"^ 

acUrif = 5) = 0.510 {Mb = 1.0 GeV ) (25) 



and 

acr^t{nf = 5) = 0.533 '^^f^ {Mb = 0.95 GeV ). (26) 

Note that the freezing value ( OBI ) turns out to be very close to that, 
phenomenologically introduced by Godfrey , Isgur in [3] to describe a 
lot of experimental data in meson sector. As shown in [17], the choice 
with Mb = 0.95 GeV = 320 MeV, and acr^t = 0.58, provides good 
agreement between experiment and BPT calculations of the splittings in 
bottomonium. Thus, in BPT the large freezing value acrit{n^f = 5) ~ 0.58 
appears to be consistent with the conventional value of ^ 230 MeV, 
a,{Mz, 2 - loop) = 0.119 ± 0.001. 

However, the large freezing value in BPT (and in QCD phenomenol- 
ogy) does not agree with the value used in the lattice parametrization 
(2) of the static potential. For example, in quenched approximation 
aiat = 0.23 was obtained in [9,10], while in BPT the minimal value acrit 
(which corresponds to the minimal value in ( |T^ . A^(min) = 218 MeV) 
is equal (y^r^{nf = 0) = 0.38, i.e., by 40% larger. Such a difference be- 
tween two numbers partly occurs due to the fact that on the lattice the 
r-dependence of the vector coupling is not seen (or neglected) at r > 0.2 
fm . 

In Table 1 the background coupling ^^(r) is compared to a^(r), cal- 
culated in PQCD, where according to the perturbative prescription the 
QCD constant A^^"* in coordinate space is defined as aJ^'^ = e^K^y^^ 
(7 is the Euler constant). These couplings (both in two-loop approxi- 
mation) appear to be close to each other only at very small distances, 
r < 0.1 fm, and have the characteristic feature that the difference 
Aa = ^^(r) — Oi^t{r) is positive at r < 0.04 fm and becoming negative 
at r > 0.04 fm; just such a change of the sign of Aa has been observed 
in lattice static potential [10]. 

Table 1: The strong coupling afj(2 — loop,r) in PQCD and 
the background coupling 0^5(2 — loop^r) in BPT 
in quenched approximation with A^^ = 237 MeV 

{A^y^ = 379 MeV, A^^"^ = 675 MeV, Mb = 1.0 GeV) 



8 



r[tm) 


p / \ 




n m 

U.Ui 




U. ioo 


n no 
u.uz 


U. iOD 


U. iOD 


n 04 

U.U4I: 


n 909 
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248 


232 


0.08 


0.301 


0.254 


0.10 


0.368 


0.272 


0.12 


0.457 


0.288 


0.14 


0.588 


0.301 



3 The static force on the lattice 

The study of the static force is especially convenient through the dimen- 
sionless function $zai(r) since this function depends on the dimensionless 
variable like ^, where a is a lattice spacing , or on ry/a (if the string 
tension is taken as the only mass scale). Recently $/at(r) was measured 
by the MILC group both in quenched case and in (2 + 1) lattice QCD 
[11] with the following results: 

First, the function <l>/at(r) = r'^Fiat{r) = {ry/a)'^F{ry/a) acquires the 
value ^latif^i) = 1-0 at the QQ separation r[: 

rjv^ = 0.769 ± 0.002 (n/ = 3) (27) 

r{v^ = 0.833 ± 0.002 (n/ = 0). (28) 

Second, the function ^lat has the value $/ai(^ov^) — 1-65 at the fol- 
lowing separation Tq (the Sommer scale) : 

rov^= 1.114 ±0.002 (n/ = 3) (29) 

ro\/^= 1.160 ±0.002 (n/ = 0). (30) 
Also in [11] the ratio 

4 = 1.449(5) (n/ = 3) (31) 
^1 

is defined with a precision accuracy. If in (27)- (31) one takes the string 
tension a = 0.20 GeV^, then the characteristic points are r{ = 0.34 fm. 
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Tq = 0.49 fm for Uf = 3 and slightly larger, r\ = 0.37 fm, Tq = 0.51 fm in 
quenched case. The numbers obtained can be used to extract the vector 
coupling ap^{r), associated with the force in lattice QCD: 

FUr)=a + ^afir), (32) 

^iatir) = ar' + ^afir). (33) 

From (27)-(31) it follows that the values of ap*{r) are equal at the points 
r{ and Vq with a good accuracy: 

afi^l^) = af{rf) = 0.307(4) (n; = 3), (34) 

afir?) = afirf) = 0.229(3) (n/ = 0). (35) 

Note that in quenched case the value (35) for ap*{r) numerically coin- 
cides with the lattice coupling aiat = 0.23 in the static potential (where 
Q^iat = const is assumed over the whole region 0.2 fm < r < 1.0 fm). 
Thus existing lattice data are consistent with the assumption that the 
derivative in this region is equal zero, a'lati'^) ~ ^■ 

The lattice number (34) in full QCD appears to be essentially smaller 
than the coupling, used in BPT and also in QCD phenomenology, in 
IR region. In order to make a conclusion, which value provides better 
description of experimental data, in Section 5 as a test we shall calculate 
the bottomonium spectrum with the lattice as well as with the BPT 
static potentials. 



4 The static force in BPT 

In BPT the static force can be presented as in (5), 

FBir)=a+^aFir), (36) 

where the coupling apir), associated with the force Fsir), can be ex- 
pressed through the known coupling asir) (10) and its derivative: 

apir) = aB{r) - ra'^ir). (37) 

Correspondingly, the dimensionless function ^^(r) is 

^B[r) = r^Fnir) = {V^rf + ^apir). (38) 
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The couphng apir) can be easily calculated from the expression ( fTUl) 
and due to negative contribution of the term with the derivative it ap- 
pears to be essentially smaller than ^^(r) (in the region r ^ 0.6 fm). 
The values of apir) for n/ = and Uf = 3 are given in Tables 2,3, from 
which the r-dependence of ap{r) is explicitly seen: 

i) At the distances r = 0.2 fm; 0.35 fm, and 0.50 fm the coupling 
ap{r){nf = 3) is smaller than ^^(r) (in the static potential) by 25%, 
18%, and 12%, respectively. 

ii) The derivative a^(r) is larger for larger Uf being approximately 
proportional /3q^ {/Sq = 11 — |n/). 

iii) At the distance r ^ 0.2 fm the coupling ap{r, nj) in BPT coincides 
with aiat both for n/ = 3 and n/ = but at larger r > 0.2 fm it manifests 
essential growth: for n/ = 3 ap{r = 0.335 fm) =0.376, ^^(^o = 0.493 
fm) =0.430, being by ~ 20% and ~ 40% larger than aiat = 0.306 from 

[11]. 

Table 2: The background couplings ap{r) = aB{r) — ra'Q{r) 
and asir) in quenched approximation {Ay^ = 379 



MeV, 


Mb = 1.0 GeV 




-- 0.419) 




r(fm) 


ap{r) 




r(fm) 


ap{r) 




0.099 


0.188 


0.272 


0.355 


0.312 


0.379 


0.118 


0.201 


0.289 


0.394 


0.324 


0.386 


0.138 


0.213 


0.301 


0.433 


0.335 


0.391 


0.158 


0.225 


0.313 


0.473 


0.345 


0.396 


0.197 


0.246 


0.333 


0.493 


0.350 


0.398 


0.236 


0.265 


0.348 


0.532 


0.358 


0.401 


0.296 


0.290 


0.365 


0.591 


0.369 


0.406 


0.335 


0.305 


0.375 


^crtit 


0.419 


0.419 



Table 3: The background couplings ap{r) = asir) — ra'Q{r) 
and aeir) for n/ = 3 (A|^^ = 370 MeV, Mb = 1.0 
GeV, a,r^t = 0.510) 
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0.490 


0.335 


0.376 


0.459 


0.591 


0.453 


0.495 



Knowing apir) the function ^^(r) (EBD can be easily calculated. Note 
that the coupling apir) (as well as Q;s(r)) weakly depends on the string 
tension through the background mass Mb = 2.236^0- (16). We have 
obtained the following numbers for the separations ri and rg, where 
<^B{ri) = 1.0 and ^^(ro) = 1.65: 

= 0.769(5) (n/ = 0) (39) 

rov^ = 1.090(5), 

and the ratio ro/ri coincides with /Tf\nf = 0) on the lattice with 2% 
accuracy: 

^ = 1.417(16)(n/ = 0). 

For n/ = 3 in BPT the characteristic separations are following, 

= 0.716(4) (40) 
rov^ = 1.044 (5) (n/ = 3) 

— = 1.458(15) 

The comparison of obtained numbers to those in lattice QCD shows 
that rl^/a{rQ^/a) in BPT is only by 8% (6%) smaller than the lattice 
values (27)-(30) while the ratio ti/tq in (40) coincides with (31) with 1% 
accuracy. 

To have a full picture - how the separations ri^/a^ tqv^ ^"^^ changing 
with a increase of flavors, below we give their values also for n/ = 5 (see 
Table 4): 

= 0.673(4), (^/ = 5) (41) 

v^ro = 1.016(5), 
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with their ratio 



1.510(16) (n/ = 5). 



Table 4: The background coupHngs aB{r) and apir) for 
n/ = 5 (A{J^ = 320 MeV, Mb = 1.0 GeV, acrit = 
0.548) 



r(fm) 






r(fm) 






0.099 


0.274 


0.381 


0.355 


0.435 


0.508 


0.118 


0.292 


0.401 


0.394 


0.450 


0.515 


0.138 


0.309 


0.417 


0.433 


0.463 


0.521 


0.158 


0.324 


0.432 


0.473 


0.474 


0.526 


0.197 


0.352 


0.455 


0.493 


0.479 


0.528 


0.236 


0.377 


0.473 


0.512 


0.484 


0.530 


0.296 


0.409 


0.494 


0.532 


0.488 


0.532 


0.335 


0.427 


0.504 


0.591 


0.501 


0.536 



Comparing (39)-(41) one can see the points ri(n/-), ro(n/) are smaller 
for larger Uf (for n/ = 5 they are by 12% smaller than in quenched 
case). Also for n/ = 5 the coupling apir) approaches a freezing value 
at smaller distances, e.g. the value apiro) /acrit is equal 0.835 (n/ = 0), 
0.843 (n/ = 3), 0.874 (n/ = 5). (See also Fig.l.) 

In contrast to lattice coupling ap*{r), which is supposed to be r- 
independent between ri and ro, (34, 35), in BPT the coupling ap{r) 
depends on r and calculated values of ap{ri), ap{ro) are given below: 



and at the Sommer scale roy/a their values are by ~ 12% larger. 



ap{rf^^/a) 
ap{rf^Va) 



0.306 (n/ 

0.366 (n/ 
0.412 (n/ 



3) 
5) 



0) 



(42) 



ap{ri^'^/^) = 0.346 (n/ = 0) 

ap{r^Q^^/^)=0A2 (n/ = 3), 
api^Q^V^) = 0.46 (n/ = 5). 



(43) 



From our analysis one can conclude that 
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Figure 1: The couplings Q;s(r) solid lines) and Q;i?(r) (dashed lines). The upper curves 
refer to case with Uf — 5 and the lower curves refer to quenched case {uf — 0). 

i) The characteristic quantities ^v^s^fo in BPT are very close to 
the corresponding lattice values being only by 8% smaller. 

ii) At the same time the coupling ai?(ro) in the "background force" 
Fb{t) is shown to be by 30%-40% larger than a^*. 

5 The 65 spectrum as a test of the vector couphng 
in IR region 

Recently it was demonstrated that the bottomonium spectrum, especially 
the splittings between low-lying levels, appear to be very sensitive both to 
the freezing value and to the r-dependence of the vector coupling [16,17]. 
Therefore just these splittings can be used for testing of Q;.sf(r) in IR 
region. For illustration we consider here three typical static potentials. 

First one imitates the lattice static potential (2) but has the form 
of the Cornell potential over the whole region r > 0. The coupling 
ck!Zai(^) — const = 0.306 is taken from lattice calculations with n/ = 3 
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from [11] which is the largest Coulomb constant obtained in lattice mea- 
surements up to now. In this potential the AF behavior of a;at(r) (which 
was observed on the lattice at small QQ separations, r < 0.2 fm [9]) is ne- 
glected and therefore such the gluon-exchange potential with aiat = const 
over the whole region is stronger than the lattice potential and gives rise 
to larger values of the splittings between levels in bottomonium. 

Second variant refers to the phenomenological Cornell potential with 
Q^st(^) = const = 0.42 which is used in analysis of charmonium [1,26]. 

Third potential is taken as in BPT for nf = 5, A^"* = 330 MeV (or 
A^(2 - loop) = 241 MeV), and Mb = 0.95 GeV. The calculations with 
first two potentials are performed with nonrelativistic kinematics while 
for third potential the solutions of relativistic Salpeter equation are used. 

The splittings between the spin-averaged masses M{nL) in bottomo- 
nium for these three potentials are given in Table 5. 

Table 5: The splittings (in MeV) between spin averaged 
masses in bottomonium for the Cornell potential 
with (A) aiat = 0.306; {B)aphen = 0.42 and for the 
BPT potential with A^^^ = 330 MeV {Mb = 0.95 
GeV) 



Splitt. 


A. 


aiat = 0.306 


B. 


C^phen = 0.42 


G. 


C^crit = 0.565 








= 0.20 GeV^ 


(7 = 


= 0.183 GeV^ 


(7 = 


= 0.178 GeV2 


exper. 




iriA 


= 4.85 GeV 


niB 


= 4.631 GeV 


rrib 


= 4.828 GeV 




2S-1S 




527 




615 




557 


563"^ 


2S-1P 




114 




97 




123 


123^^^ 


IP-IS 




413 




517 




434 


440«) 


ID-IP 




241 




260 




260 


261±2^) 


2P-1P 




359 




368 




370 


360.1±1.2 



Since the mass of r]i){nS) is unknown and in any case M{1S) < 
M(T(15')), we give here only the low limit of the experimental splitting. 

The experimental number for M(1D2) = 10161.2 zt 2.2 is obtained 
in [27]. 

From our calculations presented in Table 5 one can make important 
conclusions. First, the calculations with "the lattice" potential A with 
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o^iat = 0.306(n/ = 3) define the upper bounds of the sphttings between 
different states in bottomonium (since the AF behavior at r < 0.2 fm is 
neglected). Nevertheless even upper bounds of the 2S-1S, IP-IS splittings 
appear to be by 40-i-30 MeV smaller than the experimental numbers 
(for which we know the lower bounds since the r]b{nS) mesons are still 
unobserved, M{nS) < M{T{nS))). 

It is of a special importance to compare theoretical and experimental 
number for the ID-IP splitting which is measured now with precision 
accuracy PBj: 

A(exp) = M(1D) - M(1P) = 261.1 ± 2.2(exp) ± {th) MeV, (44) 

where M(1P) =9900.1±0.6 MeV, M(1D)= Mexp(lD2) ± (th) MeV; 

Mexp(lD2) = 10161.2 ± 1.6(exp) MeV (see a discussion of the ID -IP 
splitting in [17])- 

For "the lattice" potential A the upper hmit A (/at) turns out to be 
by 20 MeV smaller than A(exp) (the difference is about ten standard 
deviations). 

On the contrary the calculations with the phenomenological potential 
B (with aphen = 0.42) and with the BPT potential C give the precision 
agreement with A(exp) . 

In BPT the 2S-1S, IP-IS splittings are also close to the experimen- 
tal numbers being only by ~ 10 MeV smaller (the ground state mass 
M(1S)=9466 MeV is slightly larger than expected experimental number). 
The same splittings for the Cornell potential {aphen = const = 0.42) turn 
to be too large (since M(1S)= 9300 MeV is small) because the Coulomb 
part of the static potential is overestimated if the AF behavior of the vec- 
tor coupling is neglected. This AF effect is small for the ID- IP splitting 
for which agreement with experiment takes place. 

Note that in the BPT potential C the freezing value acrit = 0.565 
is essentially larger than aphen = 0.42 and to understand what kind of 
the approximation corresponds to ast{r) = const, let us introduce an 
effective coupling in BPT according to the relation: 



Our calculations of the matrix elements demonstrate (see Table 6) that 




(45) 
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(1) agff depends on the quantum numbers n, L; 

(2) the values of aeff{nL) appear to be by 30-!-15% smaUer than the 
freezing value acrit = 0.565 and those values for the IS, 2S states are 
close to aphen used in phenomenology. 



Table 6: The effective vector couphng aeff{nL) for the BPT 
potential C {K^y^ = 330 MeV, Mb = 0.95 GeV, 
a = 0.178 GeV2,aerzi = 0.565) 



state 


IS 


2S 


3S 


IP 


IP 


ID 


2D 


aeffinL) 


0.405 


0.439 


0.448 


0.495 


0.501 


0.528 


0.528 



(3) For the orbital excitations the effective coupling o^e// = 0.50 is by ~ 
20% larger than for the IS, 2S states and just this fact results in increasing 
of the splittings like 2S-1P, ID-IP which is observed in bottomonium. 

To make a decisive conclusion about the behavior of ast{r) in IR region 
it would be also important to take into account a screening of the gluon- 
exchange potential at large distances. At present there is no a theory of 
this phenomenon on the fundamental level, although in some cases the 
Coulomb screening is introduced in a phenomenological way . 

6 Conclusion 

Our study of the vector couplings ^^(r) (in the static potential) and 
apir) in the static force F{r) is performed in the framework of BPT 
with the following results. 

Due to the derivative term a^(r) an essential difference between apir) 
and asir) is observed at the QQ separations r < 0.6 fm with aF{r) 
being by 50%, 30%, 15% smaller at the points 0.2 fm, 0.3 fm, and 0.5 
fm, respectively. 

At the same time the freezing values of both couplings coincide and are 
rather large: acrit — 0.41 (n/ = 0); 0.51 (n/ = 3), and acrit = 0.58 ± 0.04 
for n/ = 5. The last number turns out to be very close to that introduced 
by Godfrey, Isgur in their phenomenological analysis. 

The dimensionless quantities rl^/o^ and ro^/a, where the function 
$(r) = r'^F{r) has the values 1.0 and 1.65, are calculated in BPT and 



17 



their values are by (6-i-8)% smaller that those calculated on the lattice 
in quenched case and in (2+1) QCD. 

In contrast to lattice observation where a^^(ri) = Q;^*(ro) = const 
and this constant is small: aiat = 0.23{nf = 0) and aiat = 0.306(n/- = 3), 
in BPT apir) at these points is found to be by 40% larger for n/ = and 
by 30% larger for = 3. Because the Coulomb constant in the lattice 
static potential is small, this potential gives essentially smaller 2S-1S, 
IP-IS, ID-IP splittings in bottomonium. 

The meaning of the Coulomb constant aphen, used in phenomenology, 
as an effective coupling in BPT is suggested. This interpretation explains 
why aphen — 0.42 is by 30% smaller than the freezing value acrit — 0.56 
for the potentials with the AF taken into account. 

The knowledge of the static force in BPT is important to perform 
precision calculations of the wave functions at the origin, hadronic decays, 
and fine structure splitting in bottomonium. 
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